Low-lying continuum structures in *B and ^Li in a microscopic model 
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We search for low- lying resonances in the *B and *Li nuclei using a microscopic cluster model and a 
variational scattering method, which is analytically continued to complex energies. After fine-tuning 
the nucleon-nucleon interaction to get the known 1^ state of *B at the right energy, we reproduce 
the known spectra of the studied nuclei. In addition, our model predicts a l""" state at 1.3 MeV 
in ®B, relative to the ^Be + p threshold, whose corresponding pair is situated right at the ^Li -I- n 
threshold in *Li. Lacking any experimental evidence for the existence of such states, it is presently 
uncertain whether these structures really exist or they are spurious resonances in our model. We 
demonstrate that the predicted state in ^B, if it exists, would have important consequences for the 
understanding of the astrophysically important ^Be(p,7)*B reaction. 

PACS number(s); 25.40.Lw, 26.65.-ht, 21.60. Gx, 27.20. -Hn 
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I. INTRODUCTION 

The ^Be(p, 7)^B and ^Li(n, 7)^Li reactions play impor- 
tant roles in nuclear astrophysics. The former process 
produces *B in the sun, which is the dominant source 
of the high-energy solar neutrinos (via the beta decay to 
^Be) H], while the latter one is a key reaction in inhomo- 
geneous big bang models |^ . In stellar environments only 
the very low-energy continuum structures of the partic- 
ipating nuclei are probed. However, terrestrial experi- 
ments often cannot study these low-energy processes di- 
rectly, and have to extrapolate the measured data (e.g., 
cross sections) from higher energies down to the astro- 
physically interesting region. Obviously, the extrapola- 
tion procedure can be strongly influenced by higher ly- 
ing structures present in the continuum. For example, 
the known l"*" and 3"*" resonances in ^B cause bumps in 
the ^Be(p,7)*^B cross section, which have to be treated 
properly if one wants to interpret the results of both the 
radiative capture Q and Coulomb dissociation [Q mea- 
surements correctly. 

The aim of the present work is to explore the contin- 
uum structure of ^B and ®Li using a microscopic model 
that has been previously applied to describe the ground 
states of these nuclei Q and in the calculation of the 
nonresonant ^Be(p, 7)®B cross section 

According to the Ref. [|j compilation, there is a res- 
onance in *B at Er — 0.63 MeV (in this paper the en- 
ergies are given in the center-of-mass frame, relative to 
the ''Be + p or ^Li -I- n threshold) and a 3^^ state at 2.18 
MeV. Many calculations predict the first state to be a 1+ 
resonance, see, e.g., Ref. [Q. In ^Li there is a bound 1"*" 
state, a 3+ resonance at 0.22 MeV, and a broad 1+ state 
at 1.18 MeV. The Ref. Q compilation does not give any 
further low-lying states in *Li and ®B. However, there are 
indications for the existence of a 2~ state in ®Li at 1.18 
MeV 0, for a 1+ state of ®B at 2.7 MeV, and for a 1" 
or 2- level in ^B at 2.9 MeV [|ll|. 

As far as theoretical models are concerned, the poten- 



tial model of Ref. seemed to indicate an additional 
Ml (1+ or 2+) state in ^B at about 1.4 MeV. However, 
that calculation was criticized |fD|l by pointing out that 
certain model assumptions in ||12[| , like the ^Be (g.s.) + p 
nature of the Ml states, were rather unphysical. Signs 
of new states have also surfaced in microscopic models. 
For example, in Ref. ||^ we noted that our microscopic 
■^He + "^He + p cluster model also seemed to indicate the 
presence of a l"*" resonance at around 1.5 MeV. 

In order to be able to make a well-founded prediction 
on the possible existence of new continuum states in ®B 
(and ^Li) in any model, one should satisfy two important 
criteria. First, the properties of the known 1+ state at 
0.63 MeV must be well reproduced by the model. In 
certain approaches, like in the potential model, this is a 
tough problem. Secondly, the recognition of any possible 
new resonance should be fairly reliable, even for broad 
states. Below we present a model which satisfies both 
conditions and predicts new states in ^B and *Li. 



II. MODEL 

We use the same model as in Refs. The eight- 

body, three-cluster C^He + '^Hc + p) wave function of ®B 
is chosen as 



E E-^{[[*?2<f'lxL(P2)],,,}, (1) 

I-J,l2 i=l 



where A is the intercluster antisymmctrizer, p2 and I2 are 
the relative coordinate and relative angular momentum 
between ^Be and p, respectively, S2 and are the spin 
of the proton and ''Be, respectively, / is the channel spin, 
and [...] denotes angular momentum coupling. While 
is a proton spin-isospin eigenstate, the antisymmetrized 
ground state {i = 1) and continuum excited distortion 
states (j > 1) of ^Be are represented by the wave func- 
tions 
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Here A is the intcrcluster antisymmetrizer between a and 
h, and are translationally invariant harmonic os- 
cillator shell model states (a = ^He, h = '^He), pi is the 
relative coordinate between a and h, h is the a — h rela- 
tive angular momentum, si is the spin of h, and Tj_^ (pi) 
is a Gaussian function with a width of jj. The Cy pa- 
rameters are determined from a variational principle for 
the ''Be energy. A similar wave function is used for ^Li 
within the '^He -I- ■^H + n cluster model space. Using (|l|) 
in the 8-nucleon Scrodinger equation, we get an equation 
for the unknown relative motions x- 

In order to ensure that all existing resonances are rec- 
ognized in the continuum, we search for the poles of 
the ^Be -I- p (^Li -I- n) scattering matrices. These ma- 
trices are generated with the help of the Kohn-Hulthcn 
variational method for scattering processes 10 . In that 
method the relative motions in are expanded in terms 
of square-integrable functions (Gaussians in our case) 
matched with the correct scattering boundary condition. 
The resulting variational scattering matrices are analyti- 
cally continued to the Riemann surface of complex ener- 
gies, using the methods of Ref. [|l5|. That is, we generate 
solutions to the scattering problem with such asymptotic 
behavior that corresponds to complex energies, 

x{e,p)^H-{kp)-S{e)H+{kp). (3) 

Here e and k are the complex energies and wave numbers 
of the relative motions, and H~ and are the incom- 
ing and outgoing Coulomb functions, respectively. These 
solutions have no physical meaning except when S is sin- 
gular, where S coincides with the physical S matrix. We 
perform this analytic continuation procedure and search 
for the poles of S. The parameters of the complex-energy 
poles are related to the resonance parameters {Er posi- 
tion and r width) through the 

e^E,~ ir/2 (4) 

equation. We mention that this way of determining reso- 
nance parameters is the most closely related to scattering 
theory [Q. The results obtained can differ from those 
coming from other methods, especially for broad states 

We note that the Be -I- p scattering picture is valid 
only below the '*He-|-'^He-|-p three-body breakup thresh- 
old, which is at 1.587 MeV experimentally (2.467 MeV 
in the case of ''Li and "'He -I- '^H -I- n). Above this energy 
three-body scattering asymptotics should be used, which 
is beyond our capacity. Although the continuum-excited 
distortion states in Eq. (|l|) can mimic three-body break- 
up at some approximate level, any states predicted above 
these energies in our model should be handled with care. 

Our approach contains a large and physically moti- 
vated part of the 8-body Hilbert space. The next step 



is to find an effective nucleon-nucleon (N-N) interaction 
that is the most suitable for the description of the prob- 
lem at hand. As we mentioned, the precise reproduction 
of the known 1+ state in ®B (and in ®Li) is a first pri- 
ority. In most of our previous calculations for ^B [Hj^] 
we used the Minnesota (MN) force. However, it turned 
out that this interaction gave an incorrect channel spin 
ratio for the 1^ state of ®Be [0, which is a member 
of a T = 1 triplet together with the l"*" states of ®B 
and ®Li. This flaw is most probably caused by the fact 
that the reproduction of this state at the right energy in 
^Be requires an exchange mixture parameter value (w) 
of the MN force which is incompatible with the experi- 
mental data in the singlet-odd N-N channel {u > 1). To 
avoid this problem in the present work, here we apply 
the modified Hasegawa-Nagata interaction jl^ , which is 
also widely used in cluster model calculations. 

We assume a pure Wigner form for the spin-orbit force 
and fine-tune the Majorana component (while adjusting 
the Wigner component to keep W + M = const.) of the 
medium range part of the central interaction, in order to 
get the 1+ state of ^B at 632 keV, its experimental value 
taken from Filippone et al. ^ (note that Ref. |^ gives 
637 keV). The required modification in the central force 
is about 5%. The strength of the short-range spin-orbit 
interaction is fixed by requiring that the experimental 
spin-orbit splitting between the 3/2^ and 1/2^ states of 
^Be be reproduced correctly. We use the same interaction 
for both *B and ^Li. Our prime target is ®B here, a 
more precise reproduction of the ^Li states and the fine- 
tuning of the N-N force for the 1+ resonance of ^Li is not 
pursued. 

III. RESULTS 

We show the spectrum of *B and *Li coming from our 
calculation in Table |, together with the experimental 
numbers. In these calculations we used = 1 inside the 
7-nucleon subsystems in (|^) in all cases, I2 = 1 ya the 
case of the 1+ and 2+ states, and Z2 = 1 and 3 in the 
case of the 3+ state. The total spin S took all possible 
values. We emphasize again that our model is not a fully 
adequate description for those states that lie above the 
three-body thresholds. These thresholds are situated at 
1.34 MeV and 2.13 MeV, respectively in our model, com- 
pared to the experimental 1.59 MeV and 2.47 MeV. One 
can see in Table ||, that the properties of the if state in 
^B are nicely reproduced. This is a significant achieve- 
ment as other models, most notably the ''Be -I- p poten- 
tial model, tend to strongly overestimate the width of 
this state. We also note that the N-N interaction which 
is fine-tuned to get the 1+ state at the right position, 
slightly overbinds the 2+ ground state. This can be ex- 
plained by the fact that the ^He, •^He)p model space is 
closer to the true wave function of the 2+ ground state 
than to that of the l"*" resonance. It means that while 
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describing the l"*" state, the N-N force has to be made 
stronger in order to compensate for those missing parts 
of the Hilbert space that cannot be represented by the 
(''He, '^He)p wave function. However, this modification of 
the force leads to the slight overbinding in the 2+ state. 

The most surprising result in Tabic || is the prediction 
of an additional 1+ state at 1.28 MeV. Note that although 
this resonance is situated rather close to our ''He+'^He+p 
three-body threshold, it is not a threshold effect or a sim- 
ilar artifact. By making the N-N force slightly stronger, 
both the energy of the state and that of ''Be is low- 
ered (while the position of the three-body threshold is 
unchanged), which makes this state more bound relative 
to the three-body threshold. The corresponding 1 J state 
in ^Li is very close to the ''Li -I- n two-body threshold, 
therefore its parameters can depend on the details of the 
model rather strongly. 

One can see in Table | that the known 3+ states appear 
in our model, although somewhat shifted to higher ener- 
gies, compared with the experimental situation. In addi- 
tion, we have a third l"*" state, although it lies above the 
three-body threshold in both *B and ®Li. Whether this 
state would appear in a realistic model, which contains 
the correct three-body asymptotics, remains a question. 

Turning our attention back to the states, we note 
that no experimental evidence has been found so far that 
would support the existence of such structures It 
seems rather unlikely, e.g., that the existence of such a 
state in ®Li, right at the ^Li -I- n threshold, could be 
reconciled with the cross section measurements of Ref. 
p9[ . Also, no indications of such states were found in 
^{t,p)^U and '^Be((i, 3He)*Li experiments Never- 
theless we feel obliged to report on the existence of these 
states in our model. Although rather unlikely, these 
states can be spurious solutions in our model. Certain 
variational approaches are known to produce such arti- 
facts |20j. Spurious states can also appear if incorrect 
boundary conditions are used. This would happen here 
also in the case of ^Li, if we searched for 2+ states in 
a three-body bound state calculation. We would find 
two such three-body bound states, one of them lying 
above the ''Li + 71 threshold. Using a wave function which 
has the correct asymptotic behavior (scattering state in 
^Li -|- n) would result in the disappearance of the sec- 
ond "state". For further details, see ^l[]. The situation 
is different in the present work. We do not know any 
conceivable reason why spurious states might appear in 
our description. Another possibility is that the predicted 
states are situated at significantly higher energies in re- 
ality, but for some reason they appear at low energies in 
our model. In the following we discuss some characteris- 
tic features of the newly found resonances, which might 
help to find them experimentally if they exist, or else 
find a way to understand them and get rid of them in 
our model if they turned out to be spurious states. 

In Table || we give the probabilities of the various an- 
gular momentum channels which build up the wave func- 
tions of the 1+ states of ®B and ^Li in several different 



forms. In these calculations a three-body bound state 
approximation was used for the 1"*' resonances, where 
the relative motions were expanded in terms of square- 
integrable basis functions. The resulting probabilities 
characterize the inner parts of the resonance wave func- 
tions rather well. 

First the weights of the orthogonal channels are 
given in Table ^ in the (L, S) and (/, ^2) represen- 
tations, respectively. Here [(^i, ^2)^, (si, S2)<S'] J and 
[[{h, si)Ij, 32]! jhjJ coupling schemes are assumed, re- 
spectively. The third way of expressing the importance 
of the various nonorthogonal channels is to calculate their 
amounts of clustering |Q . This quantity gives the proba- 
bility that the full wave function of a state lies completely 
in the given subspace (thus the sum of the probabilities 
is obviously not 1). The fourth way of characterizing the 
1+ states in Table |l| is to calculate the weights of the 
shell-model-likc I/7, j) configurations, where /y = 3/2 or 
1/2 is the spin of ''Be/^Li, while j — 3/2 or 1/2 comes 
from the coupling of the orbital momentum (Z2) and spin 
(S2) of p/n. We calculated this quantity in the follow- 
ing way. First we expressed the \I'j,pj) state in terms of 
the (L, S) components, and then evaluated its square by 
substituting the (L, S) weights. 

As one can see in Table |l|, the 1^ state receives its 
main contribution from the 1 = 2 channel spin compo- 
nents, while I2 is dominated by the / = 1 channel spin. 
In addition, as shown by the amounts of clustering, the 
1^ state receives large contributions from configurations 
which contain a ^Be in its excited state {Ir = 1/2), while 
in 1 J the role of the excited ''Be core is negligible. 

The newly found ij states, if they exist, would have 
important consequences for the radiative capture reac- 
tions ''Be(p,7)®B and ''Li(rt, 7)®Li. We do not discuss 
these consequences in detail, just mention one interest- 
ing example. In order to be able to reliably extrapolate 
the experimental ^Be(p, 7)^B cross section down to stel- 
lar energies, the precise knowledge of the nonresonant El 
cross section is necessary. The ij resonance in the vicin- 
ity of 1.3 MeV, if it exists but is not recognized, would 
change the average slope of the extracted El cross sec- 
tion. In order to show this effect, we calculated the Ml 
cross section using the above model. Due to technical 
simplifications, the asymptotic behavior of the 2+ bound 
state could only be correctly described up to 15 — 20 fm. 
Note that this is an acceptable approximation (in con- 
trast to the El cross section, where it would be totally 
unphysical) at least in the resonance regions, where the 
scattering wave functions are large at small radii. 

The resulting Ml cross section, coming from the 1+ 
initial scattering states (that is, no 3+ state is included), 
is shown in Fig. [| (solid line). Also shown in Fig. |l| is 
the El-fMl cross section (dotted curve) which we get by 
adding our Ml result to an El curve (dot-dashed line) 
which has the energy dependence of Refs. |^J^ and is 
fitted (by eye) to the low-energy Filippone data . One 
can see that, although the 1^ bump predicted by our 
model is much too strong, such a structure in this energy 
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region is not a priori ruled out by the data. A broader and 
smaller bump might be consistent with the experiments. 
Therefore, if such a state really exists, then it is weaker 
than our model prediction in Table |. We note that the 
biggest contribution to the 1^ bump in Fig. ^ comes from 
the initial scattering states where the spin of ^Be is /y — 
3/2 and the channel spin is / = 1. 

As one can observe, the average slope of the dotted 
curve in Fig. |^ is rather different from that of the El 
cross section in the 0.8 < E < 1.5 MeV region. This 
means that if the 1 J resonance really existed in then 
it would make the separation of the measured cross sec- 
tion into El and Ml components more difficult than cur- 
rently believed. 

IV. CONCLUSIONS 

In summary, we have searched for resonances in the 
and ^Li continua, using a microscopic three-cluster 
model. We used a method, the analytic continuation of 
the scattering matrix to complex energies, which is well 
suited to find any resonances, even broad structures, if 
they exist. The properties of the known 1+ states in ®B 
and ®Li are well reproduced in our model. In addition, we 
have found new 1+ resonances. We have shown that the 
state predicted at 1.3 MeV in ®B would have important 
consequences in the understanding of the ''Be(p, 7)^B re- 
action. 

As currently there is no experimental indication for 
the existence of the predicted states, we cannot be sure 
whether these structures are real or not. In order to be 
able to confirm or refute the existence of the predicted 
states, more theoretical and experimental work would be 
needed. 
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TABLE I. The low-energy spectrum of **B and *Li. The 
experimental data are taken from except for the first l""" 
state of which is from Filippone All numbers are in 
MeV. 



Theory Experiment 



r 


Er 


r 


Er 


r 


"B 2+ 


-0.215 




-0.137 




1+ 


0.632 


0.034 


0.632±0.01 


0.037±0.005 


1+ 


1.278 


0.564 






3+ 


2.98 


0.808 


2.183±0.03 


0.35±0.04 


1+ 


4.33 


1.5 






«Li 2+ 


-2.021 




-2.033 




1+ 


-0.975 




-1.052±0.001 




1+ 


0.037 


0.006 






3+ 


0.937 


0.327 


0.222±0.003 


0.033±0.006 


1+ 


2.29 


1.0 


1.18 


~1.0 



TABLE IL Characterization of the calculated 1"*" states shown in Table | by i) the weights of the 
orthogonal (L, S) components, ii) the weights of the orthogonal (1, 12) components, iii) the amounts 
of clustering of the nonorthogonal [h, /, I2) components, and iv) the weights of the shell-model-like 
{l7,j) configurations in the *B and *Li wave function. For further details, see the text. 
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IJ 




(0,1) 


0.010 


0.001 


0.791 


0.001 


0.0002 


0.737 




(1,0) 


0.188 


0.728 


0.049 


0.211 


0.726 


0.077 




(1,1) 


0.786 


0.207 


0.154 


0.759 


0.237 


0.179 




(2,1) 


0.016 


0.064 


0.006 


0.029 


0.037 


0.007 




(1,1) 


0.237 


0.888 


0.184 


0.239 


0.870 


0.142 




(2,1) 


0.465 


0.107 


0.736 


0.417 


0.115 


0.849 




(0,1) 


0.298 


0.005 


0.080 


0.344 


0.015 


0.009 




(3/2,1,1) 


0.021 


0.890 


0.113 


0.030 


0.887 


0.138 




(3/2, 2, 1) 


0.621 


0.119 


0.816 


0.611 


0.195 


0.853 




(1/2,1,1) 


0.438 


0.048 


0.050 


0.474 


0.070 


0.004 




(1/2,0,1) 


0.472 


0.017 


0.054 


0.538 


0.051 


0.008 


{h,j) 


(3/2, 3/2) 


0.108 


0.409 


0.321 


0.123 


0.406 


0.316 




(3/2, 1/2) 


0.418 


0.190 


0.317 


0.408 


0.202 


0.317 




(1/2,3/2) 


0.418 


0.190 


0.317 


0.408 


0.202 


0.317 




(1/2,1/2) 


0.049 


0.206 


0.044 


0.068 


0.191 


0.049 
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1. Astrophysical S factor for the ^Be(p, 7)**B reaction. The sohd curve shows the Ml cross 
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FIG 

section coming from the l""" initial states in the present model. The dot-dashed curve is an El cross 
section whose energy dependence is taken from [^,^ and its absolute normalization is fitted (by 
eye) to the low-energy Filippone data. The dotted curve is the sum of the two. Also shown are the 
experimental data Mm with adp = 157 mb normalization, where it applies. 
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